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Abstract 

The hole spectral function is calculated in nuclear matter to assess the relevance 
of nucleon-nucleon short range correlations. The calculation is carried out within 
the Brueckner scheme of many-body theory by using several nucleon-nucleon real¬ 
istic interactions. Results are compared with other approaches based on variational 
methods and transport theory. Discrepancies appear in the high energy region, 
which is sensitive to short range correlations, and are due to the different many- 
body treatment more than to the specific N-N interaction used. Another conclusion 
is that the momentum dependence of the G-matrix should be taken into account in 
any self consistent approach. 
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1 Introduction 


One of the main features of nuclear systems is the relevance of short range 
correlations (SRC) due to the hard core part of the nucleon-nucleon (NN) 
potential. These are relevant for the binding energy of the nuclear system in 
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general and are probed directly by the high momentum and energy region of 
the hole spectral function. In fact, SRC are entirely responsible for this region, 
where the mean field and the long range correlations can give only negligeable 
contributions. Spectral functions are extracted mainly from electron-nucleus 
scattering experiments. 

Theoretically, many approaches have been proposed for a quantitative de¬ 
scription of SRC and the corresponding spectral functions. The many-body 
aspect appears essential for a correct understanding of the correlations, and 
Bethe-Brueckner [1] or variational [2,3] methods have been extensively applied. 
Although the main features of the spectral functions in the region relevant for 
SRC seem to have a similar trend within different approaches, discrepancies 
are still present and controversial interpretations of the experimental data ex¬ 
ist in the literature. Recently [4] it was pointed out that the main features 
of the nucleon spectral function in nuclear matter can be understood in the 
framework of transport theory. In this approach the only ingredient is the 
in-medium nucleon-nucoleon cross section, which can be, in hrst approxima¬ 
tion, schematized by a single parameter throughout the entire momentum and 
energy range. The self-consistent calculation of the single particle self-energy 
and strength function then provides the nucleon spectral function. Therefore, 
phase space effects seem to dominate the behaviour of the spectral function. 

A close comparison among the results of different approaches appears nec¬ 
essary to extract the characteristics of SRC which determine the strength 
function. In particular, the dependence of the results on the particular many- 
body theory as well as on the specihc nucleon-nucleon interaction used in the 
calculations should elucidate the structure of the SRC. 

In this letter, we present calculations of the nucleon spectral function in nuclear 
matter within the Brueckner scheme using different realistic NN interactions 
and compare them with other approaches. We will concentrate the comparison 
mainly in the high energy region, where the effect of correlations is dominating. 


2 Formalism 


The hole spectral function in nuclear matter is dehned as: 
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( 1 ) 


where |4/[4) is the ground state for A nucleons with eigenvalue is 

a complete set of states for A — 1 nucleons with eigenvalues E^_i, and at is 
the annihilation operator for the normalized plane wave at momentum k. 
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In a single-particle description of nuclei only states with k < kp are occupied 
and the dual states involved in the sum are at most one hole states. In a 
correlated nucleus, the wave function of the ground state can have high mo¬ 
mentum components, and the structure of the spectral function becomes more 
complicated. 


In nuclear matter, the spectral function corresponding to the nucleon self¬ 
energy M{k,u) = V{k,u) + iW{k,u), is given by the well known result [1]: 


S{k,u) 


- lmQ{k, u) = 
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where Q{k,uj) is the single-particle Green’s function: 


G{k,uj) 
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( 3 ) 


The real and imaginary parts of the self-energy, V{k,u) and W{k,cu), are 
highly off-shell in the considered energy and momentum ranges. In nuclear 
matter, according to the BBG expansion, the whole set of two-hole line con¬ 
tributions is summed up by the diagrams depicted in Fig. 11 [1], where the 
wavy lines indicate the Brueckner G-matrix. Since it is known that correla¬ 
tions with higher number of holes give a much smaller contribution to the 
ground state energy [5-7], it is reasonable to assume that their effect on the 
spectral function is also of minor relevance, producing mainly an additional 
background. The hrst diagram (l.a) corresponds to the standard Brueckner 
approximation for the nucleon self-energy, and has the following expression 

Mi{k,u) = ^ (A;A:'|G(e(A:')+o;)|A:it')a (4) 

k'<kp 


where e{k) is the self-consistent single particle energy and the label a means 
antisymmetrization. In the energy range relevant to the hole strength function 
and for k > kp, this diagram contributes only to the real part V{k,uj) of the 
self-energy. Indeed, the imaginary part of Mi{k,Lj) comes only from on-shell 
transitions from one-particle states {k > kp) to two-particle one-hole states 
and this process is not allowed by momentum conservation. The diagram (1.6) 
takes into account transitions to one-particle two-hole states. It contributes to 
both real and imaginary parts in the considered energy and momentum range. 
Explicitely, the imaginary part W{k,uj) is given by 


W{k,u) 


^^ \{kp\G {e{h) + e{h')) 

e{p) + e{h) + e{h') - ip 
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f Ehh'p\{kp\G {e{h) + e{h')) \hh')a\^S{-u + e{p) - e{h) - e{h')) (5) 

where the sum is restricted to states h and h' with a momentum smaller than 
kp and p with momentum larger than kp. 

Once the Brneckner calculation for nnclear matter is performed, the G-matrix 
and the corresponding self-consistent single particle spectrnm are obtained. 
The G-matrix is then calculated off-shell for the entry energies needed in Eqs. 
(4,5), which provide the real and imaginary parts of the self-energy. Finally, 
the spectral fnnction is calculated from Eq. (2). 

Equation (5) is similar to the formula used in ref. [4]. The main difference 
is the inclusion in the latter of the single particle strength fnnctions in the 
phase space integral, which implies a self-consistent calculation. Fnrthermore, 
in that work the sqnare of the in-medinm scattering matrix is approximated 
by a constant average valne, which simplihes the calculation considerably. The 
average scattering matrix is then adjusted in order to reprodnce at best the 
results of ref. [ 8 ]. 

Equation (5) looks also similar to the formula for perturbative corrections in- 
trodnced in ref. [3] within the variational scheme. In this case, the g-matrix 
is replaced by the residnal interaction between hole states and two-hole one- 
particle states. The pertnrbative corrections are essential near the Fermi en¬ 
ergy since they allow the hole states to acqnire a width. However, in the high 
energy region, where we are focusing the analysis, their contribntion is only 
marginal. In the variational scheme the main contribntion comes from the 
two-hole one-particle orhtogonal correlated intermediate states, which cannot 
be cast in a formula like Equation (5). Therefore a formal comparison between 
the BBG and variational approach would be misleading. 

As a consequence the numerical comparison among these many-body theories 
appears necessary. 


3 Results 


We performed calculations of the nuclear matter spectral function with three 
different two-body potentials, the Urbana V 14 [9], the Argonne V 14 [10] and 
the Argonne vig [11]. Three-body forces were added, according to the Urbana 
IX model [12], and adjusted to reproduce the correct saturation point. We 
checked that the effect of the three-body force on the spectral fnnction is 
negligeable. The Urbana potential was used only to compare with the results 
of ref. [3] and ref. [4]. 

In Figs. 2-a and 2-b we show the spectral fnnction calculated at two values of 
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the momentum, k = 2.25 fm“^ and k = 3.5 fm“^. The three lines indicate our 
BBG calculations with different potentials as described in the figure, while 
the full circles label the results of ref. [3], where the Urbana V 14 was used. 
For both momenta the dependence on the nucleon-nucleon potential appears 
quite weak, only in some cases a discrepancy is present, which however does 
not exceed 20%. Larger deviations in the high energy region occur between 
the BBG results and those of ref. [4]. We can conclude that these discrepancies 
in the high energy region are due to the many-body treatment and not to the 
interaction employed. 

In Figs. 3-a and 3-b we compare our results (for the Urbana V 14 potential) 
to the fully self consistent calculation of Lehr et al. [4]. In this case the com¬ 
parison can be more transparent. In fact, the hrst order calculations of ref. 
[4], indicated by the dotted lines in the hgures, are exactly equivalent to our 
calculations if the G-matrix in Eq. (5) is replaced by a constant average value 
independent of momenta, directly related to in-medium nucleon-nucleon cross 
section [4]. The discrepancy with the BBG calculations is, in this case, only 
due to this approximation, i.e. to the neglect of the G-matrix momentum 
dependence. This shows the relevance of the momentum dependent correla¬ 
tions for the determination of the single particle spectral function. The fact 
that the introduction of the self-consistency moves the high energy tail to¬ 
wards the BBG behaviour appears misleading, since in the BBG calculations 
no self-consistency has been used. This observation seems equally well appli¬ 
cable to the comparison with the variational calculations [3], where also no 
self-consistency procedure is included in the many-body scheme. 

Of course these conclusions do not imply that the self-consistent procedure is 
not relevant, but only that it must be carried out with the full momentum 
dependence of the G-matrix. 


References 

[1] M. Baldo, I. Bombaci, G. Giansiracusa, U. Lombardo, G. Mahaux, R. Sartor, 
Nucl. Phys. A 545 (1992) 741. 

[2] O. Benhar, A. Fabrocini, S. Fantoni, Nucl. Phys. A 505 (1989) 267. 

[3] O. Benhar, A. Fabrocini, S. Fantoni, Nucl. Phys. A 550 (1992) 201. 

[4] J. Lehr, M. Effenberger, H. Lenske, S. Leupold, U. Mosel, Phys. Let. B 483 

(2000) 324. 

[5] B.D. Day, Phys. Rev. G 24 (1981) 1203. 

[6] M.Baldo, G. Giansiracusa, U. Lombardo, H.Q. Song Phys. Lett B 411 (1997)237. 


5 



[7] M.Baldo, G. Giansiracusa, U. Lombardo, H.Q. Song Phys. Rev. Lett. 81 (1998) 
1584. 

[8] C. Giofi degli Atti, S. Simula, L.L. Frankfurt and M.I. Strikman, Phys. Rev. C 
44 (1991) R7. 

[9] I.E. Lagaris and V.R.Pandharipande, Nucl. Phys. A 359 (1981) 331. 

[10] R.B. Wiringa, R.A. Smith and T.L. Ainsworth, Phys. Rev. 29C 1207 (1984). 

[11] R.B. Wiringa, V.G.J. Stocks and R. Schiavilla, Phys. Rev. 51C (1995). 

[12] R.B. Wiringa, V. Picks and A. Fabrocini, Phys. Rev. C38, 1010 (1988). 


6 



Figure Captions. 


Fig. 1. Diagrams contributing to the nucleon self-energy 
M{k,uj) = V{k,(jj) + iW{k,(jj). The wavy lines indicate the Brueckner G-matrix. 
(a) The first diagram is the standard Brueckner approximation Mi{k,uj) for the 
nucleon self-energy, eq.(4). This diagram contributes only to the real part V{k,Lv) 
of the self-energy, (b) This diagram takes into account transitions to one-particle 
two-hole states. It contributes to both real and imaginary parts of the self-energy 
in the considered energy and momentum range. 

Fig. 2. BBG calculations of the nuclear matter strength function S{k,iu). The solid 
line is our result with the Argonne vig potential, the long-dashed line is obtained 
with Argonne V 14 and the dot-dashed line with the Urbana vl4. The black dots are 
the variational results of Benhar et al. (a) For the momentum k = 2.25 fm“^. (b) 
For k = 3.5 fm“^. 

Fig. 3. Comparison between our results and those of Lehr et ah, ref. [4]. The dotted 
curves are obtained in that work after the first iteration step, the solid curves are 
their fully self consistent calculation; the dashed curves show the variational results 
of Benhar et al. , ref. [3]. Our BBG results calculated with the Urbana potential are 
shown here as asterisks, (a) For the momentum k = 2.26 fm“^. (b) For k = 3.59 
fm“^. 
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